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The dynamical density-matrix renormalization group technique is used to calculate spin and 
charge excitation spectra in the one-dimensional (ID) Hubbard model at quarter filling with nearest- 
neighbor t and next-nearest-neighbor t' hopping integrals. We consider a case where t (> 0) is much 
smaller than t! (> 0). We find that the spin and charge excitation spectra come from the two 
nearly independent i'-chains and are basically the same as those of the ID Hubbard (and t-J) chain 
at quarter filling. However, we find that the hopping integral t plays a crucial role in the short- 
range spin and charge correlations; i.e., the ferromagnetic spin correlations between electrons on the 
neighboring sites is enhanced and simultaneously the spin-triplet pairing correlations is induced, of 
which the consequences are clearly seen in the calculated spin and charge excitation spectra at low 
energies. 

PACS numbers: 71.10.Pm,71.10.Fd,78.30.Jw,72.15.Nj,71.30.+h,71.45.Lr 



I. INTRODUCTION 

Spin-triplet superconductivity has been one of the ma- 
jor issues in the field of condensed-matter physics. Nearly 
all the conventional and unconventional superconductors 
known to date are spin-singlet paired. The best-known 
example of triplet pairing is not a superconductor but a 
superfluid 3 He, where atomic Cooper pairs are formed in 
spin-triplet channel^ Only a few materials of spin-triplet 
superconductivity have so far been confirmed in strongly 
correlated electron systems, which iclude ruthenium ox- 
ide S^RuO^ and some heavy-fermion compounds such 
as UPt3. 3 Here, some questions will naturally arise. One 
is whether the electron correlation can take an essential 
part in superconductivity carried by spin-triplet pairs. 
Another is how the behavior differs from that of spin- 
singlet superconductivity. In this manner, research on 
spin-triplet superconductivity may offer an opportunity 
to expose unknown physical phenomena. 

Quite recently, a new mechanism of the spin-triplet 
superconductivity has been proposed in a fairly simple 
correlated electron system^ The model consistes of two 
Hubbard chains coupled with zigzag bonds and has a 
unique structure of hopping integrals: sign of the hop- 
ping integrals changes alternately along the zigzag bonds 
connecting two chains, while the sign along the one- 
dimensional (ID) chain is always negative. [A model 
where all the hopping integrals are taken to be posi- 
tive is equivalent under canonical transformation.] Under 
this sign rule of the hopping integrals, the ring-exchange 
mechanism 5 yields ferromagnetic spin correlations, and 
accordingly, attractive interaction between electrons is 
derived. In Ref. 0, the argument was developed on the 
basis of only the static properties such as pair binding en- 
ergy, spin excitation gap, and pair correlation function, 
as well as spin correlation function. Therefore, further 
investigations including dynamical properties have been 



desired. 

As for real materials, this mechanism^ may be of pos- 
sible relevance to superconductivity observed in quasi- 
1D organic conductor (TMTSF) 2 X [X=PF 6 , C10 4 , etc.], 
the so-called Bechgaard salts£i£ The system exhibits a 
rich phase diagram upon variation of the pressure and 
temperature. At low temperatures, the phase changes, 
in the order, as the spin-Peierls insulator, antiferromag- 
netic insulator, spin-density-wave (SDW) insulator, su- 
perconductivity, and paramagnetic metal, with increas- 
ing pressure. So far, experimental evidences that the 
superconducting state occurs in the triplet channel have 
been piled up, 8 although not much is known on the na- 
ture of the pairing. A newly synthesized copper-oxide 
compound P^Ba^uyOis-,^ may also be a relating sys- 
tem. This material consists of both the single CuO chains 
(as in PrBa2Cu307) and the double CuO chains (as in 
PrBa2Cu40s), and those chains are separated by insu- 
lating Cu02 plains. It has been reported that the double 
chains turn into a superconducting state below T c ~ 10 
KJ£ Although the signs of the hopping integrals seem 
not to satisfy the ferromagnetic sign rule, the structure 
of the double CuO chains bears a certain similarity to 
our modeljii 

The purpose of the present study is therefore to build 
up understanding of the ring-exchange superconducting 
mechanism by calculating dynamical quantities for the 
same model proposed in Ref. [j. To see the excita- 
tions in the spin and charge degrees of freedom sepa- 
rately, we here calculate the momentum-dependent dy- 
namical spin-spin and density-density correlation func- 
tions. We use the dynamical density-matrix renormal- 
ization group (DDMRG) method for calculating dynam- 
ical quantities?^ which is an extension of the standard 
DMRG method,^ 3 The obtained results with high res- 
olutions enable us to discuss details of the fundamen- 
tal properties on the spin and charge excitations. Thus, 
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we can find some interesting features in the low-energy 
physics of our model. 

We will show that, although the spin and charge exci- 
tation spectra are basically the same as those of the two 
weakly-coupled ID Hubbard (and t—J) chains at quar- 
ter filling, the small hopping integral between the chains 
plays a crucial role in the short-range correlations and 
low-energy excitaions; i.e., the ferromagnetic spin corre- 
lations between electrons on the neighboring sites is en- 
hanced and the spin-triplet pairing correlations between 
the electrons is induced, of which the consequences are 
clearly seen in the calculated spin and charge excitation 
spectra at low energies. We hope that the present inves- 
tigation will provide deeper insight into the mechanism 
of the spin-triplet superconductivity. 

Our paper is organized as follows. In Sec. II, we define 
the ID t-t'-U model and introduce the physical quantities 
of interest, namely, spin and charge excitation spectra. In 
Sec. Ill, exact solution of the spin (and charge) excitation 
spectrum in the noninteracting case is presented, and by 
comparing our result with the exact one, we evaluate the 
performance of our DDMRG method. We then study 
the spectra with onsite Coulomb interaction and discuss 
relevance to the spin-triplet superconductivity. We close 
with a summary in Sec. IV. 



II. MODEL AND METHOD 




FIG. 1: Schematic representations of (a) the lattice structure 
of our model and (b) noninteracting band dispersion at t' /t > 
1. 

We consider the ID Hubbard model defined by the 
Hamiltonian 

i,<r 

+ i'E(c] +2CT c iCT + H.c.) + E/En iT n a . (1) 



where c\ a (ci a ) is the creation (anihilation) operator of 
an electron with spin a (er =|, \) at site i, ni a = c| CT Ci CT 
is the number operator, t' and t are the nearest-neighbor 
and next-nearest-neighbor hopping integrals respectively, 
and U is the on-site Coulomb interaction. We choose the 
signs of the hopping integrals t and t' to be positive, 
so that the two spins on a triangle of the lattice align 
ferromagnetically when U is large.— The lattice structure 
of our model is shown schematically in Fig.QJa). We call 
the chain along the t (t 1 ) hopping integral the i-chain 
(i'-chain). The dispersion relation is given by 



£fe = 2t cos ka + 2t' cos 2ka, 



(2) 



where a is the lattice constant along the i-chain (we set 
a = 1 hereafter). We consider the case where there are 
four Fermi momenta ±fci?i and ±kp2 (|^F2| > |); i.e., 
the case where there are two branches in the noninter- 
acting band dispersion [see Fig. [l{b)]. This case occurs 
when the hopping integrals satisfy the condition 



i> cos 2 [(2 - p)tt/2] 
t sin' 



(3) 



where p is the band filling. In this paper, we restrict 
ourselves to the case where t' is a few times as large as 
t and the system is at quarter filling, p = 1/2. Hence, 
the model can be regarded as a double i'-chain Hubbard 
model weakly coupled by the t-chain. 

Because t is much smaller than i! , it would be very use- 
ful to allow a case of t — for familiarization with our 
results. In the limit of t — ► 0, the system is equivalent to 
the two independent ID Hubbard chains at quarter filling 
since electrons are distributed equally to the two chains. 
The noninteracting band dispersion reads Ek = 2t' cos 2k 
and there are four Fermi momenta ±fcpi = ±57r/8 and 
±&F2 = ±37r/8 in the original Brillouin zone defined for 
the t-chain. We define 2kp = fcp2 — &fi = tt/4 as the 
'nesting' vector in our model at t — 0. We use this defi- 
nition of the Brillouin zone throughout the paper. 

We calculate the spin excitation spectrum 
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defined with s+ = (1/vT) £ r e lqr c\ |C r |, and the charge 
excitation spectrum 
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defined with n q — {^/VL)^ rij e lqr n r( j. Here, \^>q) and 
Eq are, respectively, the wavefunction and energy of the 
ground state of the Hamiltonian Eq. U). The DDMRG 
technique is applied to calculate the excitation spectra. 
We here use the open-end boundary condition (OBC) 
for accurate calculation, because the system is relatively 
hard to deal with by the DMRG method due to large 
long-range hopping integrals. 13 When the OBC is used, 




FIG. 2: Exact spin excitation spectrum 2S(q, uj) [= N(q, uj)] 
for t = 0.25 and t' = 1 in the noninteracting case (U — 0). 
Broadening of the spectrum r\ = 0.03 is introduced. Five 
momenta indicated with arrows in the left side denote gapless 
points. Inset: Exact spin structure factor S(q). 



FIG. 3: Spin excitation spectrum 2S(q,u>) [= N(q, u>)] in 
the noninteracting case (U = 0) calculated by the DDMRG 
method. The same set of the parameter values as in Fig. [2] 
is used. Right (left) panel shows the result for L — 24 and 
r\ — 0.1 (L — 48 and r\ — 0.04). Inset: Spin structure factor 
S(q) obtained from cj-integration of S(q,ia). 



we need to use the quasimomcnta k = nm/(L + 1) for 
integers 1 < m < L on a chain with L sites in order 
to express the momentum-dependent operators s+ and 



III. RESULTS 

Noninteracting spectrum 

First, let us consider the noninteracting case, U = 0, 
where the model is exactly solvable. In this case, an 
excitation corresponds to creating a particle-hole pair for 
the ground state and therefore we can obtain the exact 
spectrum of spin excitations: 

S(q,uj) = lim — — 

(6) 



where small r\ is introduced to regularize the poles at par- 
ticular frequencies ui. Note that the spectrum of charge 
excitations at U = is exactly twice as large as that 
of the spin excitations, i.e., N(q,u>) = 2S(q,w), because 
N(g, uj) is just a sum over both up and down spins. 

In Fig. [2l we show the exact noninteracting spin ex- 
citation spectrum S(q,u>) given in Eq. ([6]). For small 
t/t' (< 1), the spectrum contains two predominant fea- 
tures: (i) large-weighted structure consisting of the dou- 
ble sine curve, whose dispersions are approimately writ- 
ten as uj ~ (At' ± 2t)sinq, and (ii) small- weighted con- 
tinuum structure at low frequencies, which arises from 
excitations between different branches of the noninter- 
acting bands. A zero-energy excitation is caused by the 
creation of a particle-hole pair just at the Fermi level sp, 
so that the gap closes at five momenta q = 0, k^2 — &fi 
(= 7r/4), 2/cpi, 2/cpi, and fcpi + fcp2- 

Now, using the DDMRG method, we attempt to repro- 
duce the noninteracting spin excitation spectrum. The 
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result is shown in Fig. [2] Since the noninteracting model 
poses a nontrivial problem to the DDMRG technique, 
it gives us a relevant accuracy test. When carrying out 
the DDMRG calculation, we have to take into account 
the required CPU time. Usually, the DDMRG method 
takes much longer CPU time than the standard DMRG 
method because the excited states must be obtained and 
an asked quantity must be calculated (almost) individ- 
ually for each frequency. Additionally, a required CPU 
time tcpu increases rapidly with frequency to and system 
size L in the DDMRG calculation, which is estimated ap- 
proximately as tcpu «u° (1 < a < 2) and as tcpu oc L 
if we keep other conditions. Hence, it would be efficient to 
take a relatively small system for obtaining an overview 
of the spectrum and a larger system for studying the de- 
tailed structure at low frequencies. 

Let us then check our DDMRG result with the exact 
spectrum. In the right panel of Fig. [3l we show the spin 
excitation spectrum S(q,u>) at U = calculated with 
the DDMRG method in a chain with L = 24 sites. The 
structure of the double sine curve can be clearly seen. 
However, it is hard to see dispersive structures in the 
continuum spectra at low frequenciese because only dis- 
crete peaks can be obtained due to the finite-size effect. 
We need to take larger systems to resolve this problem 
since the resolution of spectrum can be improved in pro- 
portion to the system size. We therefore choose to double 
the system size, L = 48, and consider the low-energy ex- 
citations. The result is shown in the left panel of Fig. [3] 
The resolution is obviously improved and we can now 
confirm the five momenta at which the zero-energy ex- 
citations occur. Moreover, we can see good agreement 
in the spin structure factor S(q) [= uj )] between 

the DDMRG and exact results, as shown in the inset of 
Fig. [3J Thus, we are confident that our DDMRG method 
indeed enable us to study the detailed structures of the 
relatively complicated spectrum. 

For the information, we keep m = 400 (800) density- 
matrix eigenstates to obtain the spectrum for L = 24 (48) 
sites. Note that a larger m value should be necessary to 
get the true ground state and excited state sZ q \^o) ( or 
n_g|^o)) °f the system. We therefore set m = 1200 in 
the first 4-5 DDMRG sweeps. 



Spin excitation spectrum 

Next, let us see how the spin excitation spectrum is 
modified by the inclusion of the onsite Coulomb interac- 
tion U. In Fig. H we show our DDMRG result for the 
spin excitation spectrum S(q,ui) at t — 0.25, t' — I, and 
U = 10, where we use the clusters with L = 24 (right 
panel) and L = 48 (left panel). 

Roughly speaking, the lower edge of the spectrum con- 
sists of three sine curves with four nodes at q ~ 0, 
^F2 — kpi (= 7r/4), 7r — (fcp2 — fcpi) (= 37r/4), and w. The 
excitation gap seems to close around these nodes. This 
is consistent with previous theoretical studies^-i 5 which 




FIG. 4: Spin excitation spectrum S(q,ui) at t = 0.25, t' = 1, 
and U = 10 calculated by the DDMRG method. We use the 
clusters L = 24 with the broadening rj = 0.1 (right panel) 
and L = 48 with rj — 0.04 (left panel). Inset: Spin structure 
factor S(q) obtained from tj-integration of S(q,uj). 



have suggested that no spin gap exists in the strong cou- 
pling regime. On the other hand, the higher edge of the 
spectrum is approximately represented as a sine curve 
lu ~ sing as in the case of U = 0, which comes from 
the creation of the particle-hole pairs within the same 
branch. 

Let us now take a closer look at the spectrm. As far 
as the spin degrees of freedom are concerned, the model 
Eq. {T]) for large U may be mapped onto a two-chain 
t— J model coupled with zigzag bonds. For small t/t', 
the antiferromagnetic interaction along the £'-chain, J', 
must be much larger than that along the i-chain, J, if 
we assume that the exchange interaction comes from the 
second-order perturbation of the hopping integrals with 
a fixed U, i.e., J' (~ t' 2 ) > J (~ i 2 ). Consequently, 
the features of the spectrum can be basically interpreted 
as those of the ID quarter-filled t— J modeL 17 ' 18 ! 19 The 
nodes at the lower edge of the DDMRG spectrum occur 
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at the momenta q — 0, 2k F , 2tt — 2fc F , and tt. Also, we 
can see considerable enhancement of spectral intensities 
around q — tt/4 (— fcp2 — ^fi) in comparison with the 
nonintcracting spectrum. This indicates that the onsite 
Coulomb interaction enhances the antiferromagnetic cor- 
relation with a period of 4 times the lattice constant along 
the i'-chain, which can be easily expected from the fact 
that the 2fc F -SDW correlation is the most dominant for 
small J in the ID t— J model at quarter filling. This re- 
sult may be compared with the 2/cp-SDW state observed 
experimentally in (TMTSF) 2 X.i& 

We then examine the effects of a small hopping in- 
tegral t, which leads to the antiferromagnetic interac- 
tion J along the t-chain as mentioned above. From the 
viewpoint of the spin degrees of freedom, magnetic frus- 
tration must be brought because triangular lattices are 
formed of only the antiferromagnetic interactions. Al- 
though J is much smaller than J', we can clearly see 
the influence in our DDMRG spectrum; there are two 
nodes around q = 37r/4, i.e., q ~ 2kpi (> 37r/4) and 
2tt — 2fcp2 (< 37r/4). These nodes are collected into a 
single node at q = 37r/4 when t = 0. This split actually 
signifies a tendency to a formation of the 2fcp-SDW state 
along the i-chain as well as to a collapse of the 2fcp-SDW 
state along the i'-chain. With increasing t/t', the node 
at q = 2fcp2 approaches q = tt/2 and the adjacent spec- 
tral weight increases, whereas the node at q — 2fcpi goes 
away from q = 3tt/4 and the weight goes to zero. In 
other wards, the hopping term t weakens the 2fcp-SDW 
oscillation along the i'-chain since the competing antifer- 
romagnetic correlation along the i-chain is enhanced. In 
fact, the spectral weight around q = 7r/4 will certainly 
diminish as t/t' increases. 

Another noticeable feature is that the spectral weight 
around q = Sir /A is obviously smaller than that around 
q = 7r/4, while the spectrum should be symmetrical 
about q — tt/2 in the case of t — 0. This result im- 
plies the presence of ferromagetic spin fluctuations at low 
energies. For clearer understanding, we study the spin 
structure factor S(q). As shown in the inset of Fig. [H 
it is evident that S(q) around q — 7r/4 is greater than 
that around q = 3tt/A; otherwise, S(q) should be almost 
symmetrical about q = tt/2. This result indicates an en- 
hancement of the ferromagnetic spin correlation between 
the neighbouring sites along the i-chain. We may explain 
this as follows. 

A real-space behavior of the spin-spin correlation may 
be derived from the Fourier transform of S(q), 

(^> = i^%)e^-^), (7) 
q 

with Sf — (riii—nii)/2. We find that the spin correlation 
along the i'-chain is not affected so much by a small t 
value but the decay length of the 2/cp-SDW oscillation is 
slightly shortened. On the contrary, the t hopping term 
plays a prominent role in spin correlation between the 



i'-chains, which may be estimated as 

(S*S* +R ) cx cos [(7r/4)i?] x (decaying term), (8) 

where R is an odd number. For R = 1, we obtain 
(SfSf +1 ) = 0.00745, which indicates the presence of 
ferromagnetic correlation between two electrons at the 
neighboring sites. This result is consistent with a sce- 
nario of spin-triplet superconductivity where the pairing 
of electrons occurs between the inter i'-chain nearest- 
neighbor sites, which is proposed in Ref. Q. 

Charge excitation spectrum 

Finally, we study the charge correlation function in 
the presence of the onsite Coulomb interaction U. In 
Fig.[5j we show the charge excitation spectrum N(q, uS) at 
t = 0.25, t' = 1, and U = 10, calculated by the DDMRG 
method for the chains with L = 24 (right panel) and with 
L = 48 (left panel). The outline of the spectrum, which is 
roughly expressed as uj w 4t' sin q, seems to be similar to 
that in the noninteracting spectrum. This result reflects 
the fact that the overall dispersion is hardly affected by 
the onsite Coulomb interaction. However, we notice that 
two distinct features emerge in the low-frequency range 
of the spectrum, which are discussed below. 

One is the increase in the low-energy spectral weight 
around q — tt/2 in the excitation spectrum. The exci- 
tation gap seems to close there, so that the lower edge 
of the spectrum consistes of the two sine curves with 
three nodes at q ~ 0, tt/2, and tt. We may thus have 
one gapless charge mode. This result reflects a tendency 
toward the Peierls instability, i.e., the formation of the 
4/cF-charge-density-wave (CDW) along the i'-chain. We 
find that the spectral weight for the 2fcp-CDW correla- 
tion is instead much reduced. This feature would be more 
evident if we observe the momentum-dependent charge 
structure factor N(q), which is shown in the inset of 
Fig. [51 We find that N(q) takes the maximum value 
at q = tt/2 and goes down along the practically straight 
lines to q = and tt. This dependence is consistent with 
results in the ID quarter-filled Hubbard chain at large 
U. 20 ; in fact, the result is almost equivalent to the result 
for noninteracting spinless fermions at half filling, as far 
as the charge degrees of freedom are concerned. 

The other is the apperance of large-weighted sharp 
peaks around q — at u) 0. The point is as follows; the 
spectral weight of the peaks around q = is larger than 
that around q = tt in the low-frequency range, and they 
are also gathered at lower frequencies. This implies that 
the electrons tend to come in the neighboring sites along 
the t-chain, which is associated with the pairing of two 
electrons between the inter t'-chain sites; accordingly, the 
pairs tend to be arranged alternately along the t'-chain. 
Note that the two momenta q = and tt should be equiv- 
alent when t — 0, so that the present result implies that 
the small hopping integral t enhances the pairing fluctu- 
ations at low energies. 
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3tt/4 = = 
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FIG. 5: Charge excitation spectrum N(q, ui) calculated by the 
DDMRG method for the same set of the parameter values as 
in Fig. [4] We use the clusters L — 24 with broadening rj = 0.1 
(right panel) and L = 48 with rj = 0.04 (left panel). Inset: 
Charge structure factor N(q) obtained from w-integration of 
N(q,co). 



Additionally, we can estimate the so-called Tomonaga- 
Luttinger liquid parameter K p from the derivative of 
N(q) at q = 0, 



ir dNjq) 



(9) 



whereby we find the value K p ss 0.637. Because, in the 
presence of one gapless charge mode, the pairing cor- 
relation is dominant in comparison with the 4Ap-CDW 
correlation when K p > 0.5 ) 21 ' 22 our estimated value of 
K p is consistent with the occurence of spin-triplet super- 
conductivity proposed in Ref. |4j. 



IV. SUMMARY 

We have calculated the spin and charge excitation 
spectra of the two-chain zigzag-bond Hubbard model at 
quarter filling in order to seek for consequences of the 
spin-triplet pairing induced by the ring-exchange mech- 
anism. The model is topologically equivalent to the ID 
Hubbard model with nearest-neighbor t (> 0) and next- 
nearest- neighbor t' (> 0) hopping integrals. We here 
have considered the case at t C t'. We have used the 
DDMRG technique to calculate the excitation spectra. 

We have first demonstrated the accuracy of the 
DDMRG method by reproducing the noninteracting ex- 
act spectrum. We have suggested that, for practical 
calculations, it is neccesary to adopt a relatively small 
system for obtaining an overview of the spectrum and 
a larger system for investigating detailed structures of 
the spectrum at low frequencies because the required 
CPU time increases rapidly with increasing the frequency 
and/or system size. Thus, the DDMRG method enables 
us to study the details of relatively complicated struc- 
tures of the spectra. 

Then, we have investigated how the spectra are de- 
formed when the strong onsite Coulomb interaction U 
sets in. We find that the spin and charge excitation spec- 
tra are basically the same as those of the ID Hubbard 
(and t-J) chain at quarter filling; i.e., the spectra come 
from the two nearly independent ID chains where the 
2fcp-SDW and 4fcp-CDW correlations along the i'-chains 
are enhanced. However, we find that the hopping inte- 
gral t plays a crucial role in the short-range spin and 
charge correlations; i.e., the ferromagnetic spin correla- 
tions between electrons on the neighboring sites is en- 
hanced and the spin-triplet pairing correlations between 
the electrons is induced, of which the consequences are 
clearly seen in the calculated spin and charge excitation 
spectra at low energies. Our DDMRG calculations for the 
spin and charge excitation spectra have thus supported 
the ring-exchange mechanism for spin-triplet supercon- 
ductivity where the pairing of electrons occurs between 
the two chains. 
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